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Abstract 

The electrical resistivity, thermoelectric power and electronic thermal con- 
ductivity of simple (isotropic) metals are studied in a uniform way. Starting 
from results of a variational solution of the Boltzmann equation, a gen- 
eralized Matthiessen rule is used in order to superpose the inelastic (or 
not) electron-phonon and elastic electron- impurity scattering cross sections 
("matrix elements"). The temperature dependence relative to these pro- 
cesses is given through simple functions and physical parameters over the 
usually investigated range of temperature for each transport coefficient. The 
coherence of such results is emphasized. 

PACS05.60.+W, 17.15. Jf, 44.90,+c 



I. INTRODUCTION 

To consider the most simple models is sometimes necessary for studying the charge 
carrier and thermal transport in solids. Indeed mathematical approximations used on 
more sophisticated models may make less evident, or even completely suppress, certain 
effects. For the transport coefficients the latter statement surely concerns the thermoelec- 
tric power and the thermal conductivity. For studying the electrical resistivity theoretical 
methods are surely simpler, better developed and more reliable. 

When there is more than a single scattering source the electrical resistivity is suf- 
ficiently accurately treated in metals under the assumption of Matthiessen ruleEHl (see 
exceptions in Concluding Remarks). Further assumptions are necessary to find simple 
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and effective expressions for the thermoelectric power and the thermal conductivity. The 
Gorter-Nordheim rule!, which is usually applied to treat the thermoelectric power, is 
based on the assumption of the Mott (-Jones) formula validity, however, restricted to the 
elastic scattering! and the adiabatic approximation!!!!. For the thermal conductivity the 
Wiedemann-Franz law is a usual basis on which the separation of the electronic and 
lattice contribution to the experimental thermal conductivity is made, whence for com- 
paring experimental and theoretical data. The validity of the Wiedemann-Franz law is 
also restricted to elastic collisions. Clearly to resolve the question of elasticity and inelas- 
ticity of the scattering is essential for comparing experimental data and theoretical laws 
in particular for the thermoelectric power and the thermal conductivity. 

The most simple inelastic scattering model of the charge carriers by phonons describes 
most features! of the thermoelectric power in normal metals. Consequently, an approach 
well accounting for the energy transfer between the charge carriers and the scattering 
quasiparticles should be applied beyond that based on the Mott formula! . 

It is known that the elastic approximation is sufficient for considering the scattering of 
the conduction electrons by static impurities. This leads to a constant, temperature inde- 
pendent residual resistivity, when the potential is appropriately chosen and the electron 
gas is degenerate. The simplest and well founded potential which leads to the residual 
resistivity in metals is the screened Coulomb potential (or Yukawa potential) provided 
that impurity concentration is sufficiently low. 

Phenomenological consider at ionsEttj show that without assuming a behavior for the 
electrical residual resistivity one cannot explain the low temperature behavior of the elec- 
tronic contribution to the thermal conductivity. At moderate temperatures, along with the 
scattering by phonons the impurity scattering explains the magnitude and the basic fea- 
tures of the temperature dependence of the thermal conductivity of normal nonmagnetic^ 
and even magnetic0 metals, provided that the magnetic moment in the latter is periodic 
and well localized. It was shown indeed in Ref. devoted solely to the thermal conduc- 
tivity, k, that the magnetic scattering was of equal importance to the impurity and the 
phonon scattering. However, k was neither compared to the electrical resistivity p, nor to 
the thermoelectric power S, for example calculated in the same temperature interval in 



Ref. [12]. Moreover the final results were not presented through standard functions such as 
the Bloch-Griineisen ones. We shall do so in the present paper, but confining ourselves 
to non magnetic impurity and phonon scattering only. 

As in our previous papers we shall apply the Ziman variational approach0 within the 
Boltzmann equation formalism, an approach in which the notion of "relaxation time" 
is not used. Thus, the relationship to the Mott formula is not straightforward!. How- 
ever, terms responsible for elastic and inelastic scattering can be easily selected, at least 
when the conduction band is parabolic. Within this approach one can generalize the 
Matthiessen rule: we consider a superposition of the transition probabilities themselves 
within the Boltzmann equation formalism!lll instead of the (integrated) relaxation time 
components. Therefore the electrical resistivity, the thermoelectric power and the thermal 
conductivity are considered here with respect to impurity and phonon scattering on the 
same ground, i.e. both model and approximations. For this purpose it is sufficient to 
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introduce a single parameter which accounts for the mutual contribution of the phonon 
and the impurity scattering to the transport coefficients. 

A similar model accounting for the scattering by phonons and impurities was consid- 
ered in the fiftiesS within the Kohler variational formalism00. However, the mathematical 
method was more complicated and less physical than that developed later by ZimarH 
It did not allowHS to get qualitatively correct results for the thermoelectric power in 
contrast to those for the electrical resistivity and the electronic thermal conductivity. 

The finding of an accurate solution of the (Bloch-)Boltzmann equation is not a past 
problem of the fifties or the seventies but still a present-day problem. Numerical methods 
are used for this purpose. Recently^ the Allen methodlll of solving the equation with a 
good numerical convergence has been used to study the influence of the electron-electron 
scattering on the electrical and thermal conductivity. A comparative discussion of the 
approximation made while using the Ziman variation methodS and the Allen oneS is 
given in Sec.V. 

In the following Section II we shortly describe the Ziman variational method and the 
general results when applied to the model here defined. Certain general results and for- 
mulae are also given in Appendix A for completeness. In Section III the matrix elements 
describing the phonon and impurity scattering are calculated and final formulas for the 
electrical resistivity, electronic thermal conductivity and thermoelectric power are system- 
atically presented. In Section IV the temperature dependence of the transport coefficients 
and that of the Lorenz number is numerically illustrated. The visual presentation is made 
by using dimensionless functions. An analysis of appropriate material constants allow- 
ing to estimate the magnitude of the final transport coefficients is found in Appendix B. 
Section V contains the concluding remarks. 



II. MODEL, METHOD AND APPROXIMATIONS. 

We assume that the current carriers are free electrons with a parabolic isotropic energy 
band e(k) = fr 2 k/2m spectrum. We also assume that the concentration of the electron 
gas (of both subbands with spin up or down) n e = (l/3n 2 )(2m/h 2 ) 3 / 2 e F 3 / 2 (or the Fermi 
energy e F ) is sufficiently high to consider the gas to be degenerate. Then one can assume 
that the chemical potential is C(T) — £ v, an d all relevant functions of z — (/k B T can be 
considered as asymptotic expressions of the argument e F /k B T . In particular the Fermi- 
Dirac integrals 

F n (z) = [°°dx — — with z = (1) 

Jo 1 + exp[x — z\ k B l 

can be represented as 

F n (z) = z n+1 /(n + 1), with z — £ F j k B T (2) 
for sufficiently large z. 
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The carriers are assumed to be scattered by static, nonmagnetic impurities, which 
concentration is sufficiently low to have a linear influence on the electrical resistivity . The 
scattering of carriers by phonons will be considered on the same footing. However, the 
phonon (quasiparticle) scattering will be treated as inelastic while the impurity scattering 
is elastic in its nature, since the impurity is usually assumed to be more heavy than 
the electron effective mass. Even though some calculations can be performed without 
assuming the form of the potential) we assume the form of the impurity potential to be 
the screened (Yukawa) potential 

U(r) = _^ 2 exp(-Ar) (3) 
r 

with r = |r|, Z — Z impurity — Z hostion the effective valence and the screening constant 
A. The Fourier transform U(q) = 4^Ze 2 /(X 2 + q 2 ) of the above function is responsible 
for the magnitude of the carrier scattering; q with q = |q| is the electron wave vector 
transferred during the scattering. The transport relaxation time following from U(q) - 
and consequently the residual resistivity - is temperature independent for a degenerate 
electron gas. Their magnitude is related to the Fermi wave vector k F or the energy e F 
through the function! 

2 C (A 2 d +g 2)2 = ln (! + 4 *4) - 4*4/(1 + 4*4) = <P(2k FX ) } (4) 
where k FX = k F /X . 

Note that the standard consider at ionsa yield that A 2 ~ n e 1//3 and thus one can estimate 
2k F = (3n) 1 ^ 6 (aHk F ) 1 ^ 2 = 1 since the effective Bohr radius an = h/(me 2 ) is of the order 
of k F . 

The phonon system and the electron-phonon scattering are considered in the Bloch 
approximation, i.e. phonons in the harmonic approximation and the electron-phonon 
interaction in the deformation potential approximation. We restrict the magnitude of the 
phonon wave vector, q < q^, to be within the Debye sphere. Then, for sufficiently high 
values of the Fermi energy, z =— > oo (see, e.g., Appendix B in Ref. |8|), the phonon limited 
transport coefficients are expressed in terms of the Bloch(-Gruneisen) functions 

i//t 

*">=/ («-l)(l-e~) - <5) 

A relation between the deformation potential approximation and the method based 
on the dynamical structure factor introduced by van Hove§ is described in Sec. V. 

As we mentioned in the Introduction we apply the Ziman variation method^ to con- 
sider the combined effect of the impurity and phonon scattering on the transport co- 
efficients. The method consists in calculating the elements of a matrix describing the 
scattering and two vectors governing the flow of the electric charge and heat. The "ex- 
act" results correspond to an infinite order matrix and infinite dimension current vectors 
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(see, e.g. Appendix A and Sec. IV in Ref. |5[). Here, we truncate the matrix and vectors to 
four elements and two components, respectively. This is the lowest order approximation, 
which accounts well for both the electrical and thermal current if the variational trial 
functions are appropriately chosen00. This choice is known to be 

^(k) = (va)( £ -Cr\ (6) 

where % — 1 and 2, while a denotes the unit vector in the direction of the applied 
temperature gradient or that of the external electric field causing the charge flow, and 
v = {Ti/m)Vk£ is the electron velocity. 

In this approximation the transport coefficients : the electrical resistivity p, the ther- 
moelectric power or the Seebeck coefficient 5* and the (electronic) thermal conductivity k 
for cubic crystals, are given as such 

p = (p 11 p Z2 -p* t )n-\ (7) 

S = [P 22 JiU 2 - P 12 (JiU 2 + J 2 U ± ) + P U J 2 U 2 ] (Til)- 1 (8) 

K = (Jlul + J%Ul - 2J 1 J 2 U 1 U 2 ) (Til)- 1 . (9) 

with 

n = p 22 ji - 2 p 12 j ± j 2 + p ±1 r 2 , (io) 

the elements of the scattering matrix being 

p l3 = jdk ydk'c^k'K^k'), 

Uij (k, k') = [^(k) - ^(k')] [^-(k) - ^(k')] , (11) 
and the components of the trial currents [] 



J l = -e Jdk (-^j (va)A(k) , 



(12) 



Ui= /dk (-^\ (va)(e-C)^(k) , (13) 



1 The factor (v. a) in the expression for the trial currents in Ref. 5 was inadvertently omitted. 
In the trial functions (12)-(13) this factor is used as in Ref. 2, and in our previous papers where 
it was denoted (k.a). 
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where C (k, k') is the transition probability to be determined on the ground of an appro- 
priate Hamiltonian describing the scattering. The function / denotes the Fermi-Dirac 
one, / = 1/[1 + exp(e — ()/k B T], and — e(e > 0) is the electron charge. 

The integrals in (|TT|) and fll^)-fll3D should be taken over the Brillouin zone but in our 
free electron considerations we let 

/ dk ...= 2 7 /-/dS,\ 

J (277)3 7 J\\V k e\) y } 

where dS^/jV^I is an element of the isoenergy surface. The integrals with respect to the 
electron energy in (10) for the degenerate electron gas are described in Appendix B of 
Ref. 8. In the case of elastic scattering by impurities, as here, the integrals can be treated 
as simply. Namely, we take advantage of the approximate equalities 

oo 

J d£(-d//dfi)(e - 09(e) = (7f 2 /3)(k B T) 2 [dg/de} c , (15) 
o 



J de(-df/de)(e - () 2 g(s) = (vr 2 /3)(k B T) 2 g(() . (16) 
o 

where g(e) is an arbitrary function. 

The trial currents (|I~2|)-(|1~3|) can be represented! as acombination of the Fermi- 
Dirac integrals (see Appendix B) : J x = -J (k B T) 3 / 2 L (z), J 2 = - J^Tf' 2 L^z) , 
Ui = -UoikvTf^L^z), U 2 = -U {k B T) 7 / 2 L 2 (z), where U and J = eU are material 
constants. 

For the electron-phonon system the P-matrix elements can be derived (see Ref. |8]) 
without any further approximation. They are expressed by double integrals with respect 
to the electron energy and the transferred electron momentum (wave vector). For a Debye 
cutoff and in the limit z — » oo they can be represented^ by P^- = P (k B Ty + i~ 2 Vij, where 
Po is a material constant and Vij = Vij(t) are dimensionless combinations of the the 
Bloch-Griineisen functions of the argument t = T/Tp (Tu is the Debye temperature). 

For the impurity scattering an expression similar to the above is readily obtained. In 
this case the V%j 's are simply related to <P(2k F /\) defined in (4). The explicit form of the 
functions Vij will be given in the next Section. The magnitude of Po and that of other 
material constants are estimated in Appendix B. 

In view of the above, by calculating the scattering matrix elements in the degenerate 
limit, retaining only V\\(^> P12/P22) in case of the electrical resistivity and the lowest 
order terms in the 1/z in the trial currents (12)-(13), we obtain 



2 For further convenience, the factor t 5 used in Ref.8 is introduced here directly in the definition 



of V^. 



Pii PoPn 

72 J2 F 3 ' 



S = S (t) + ^(f) — 

Bp 



1 _ P22 _ Pq TV22 

k ~ {^/?,)C T{k^TfJl ~ U$ (7r 2 /3) 2 (A; B T)% F 3 



(17) 



(19) 



where J\ in ([[7]) and ([T^) is the trial current in its asymptotic form and C = (n 2 /3)(k B /e) 2 
is the Sommerfeld value of the Lorenz number. The explicit form of the functions So(t) 
and Si(t) of (pT8|) - having a more complicated dependence on the matrix elements Vij - 
are discussed in Section III B. 



III. SCATTERING MATRIX CONTRIBUTIONS. 

As seen from (|TT|)-(|T9|) it remains to calculate the elements of the P-matrix in order 
to find the final expressions for the transport coefficients. In our case the crux is to take 
into account the combined effect of the scattering by impurities and phonons. 

As we mentioned in the Introduction we apply the Matthiessen rule in a more general- 
ized wayBOt^ i.e. to the microscopic transition probabilities C(k,k') themselves appear- 
ing in Eq.([U]) rather than to the final coefficients: the electrical and thermal resistivity. 
The variational method enables us to use the same approximation to both these coeffi- 
cients and the thermoelectric power through the integrated probabilities, i.e, we assume 

C(k, k') = C {mp) (k, k') + C (ph) (k, k') . (20) 

This means that the resultant scattering probability for the phonon and impurity scatter- 
ing can be simply considered as a sum of the scattering probabilities of the impurities (the 
phonon scattering being quenched) and by phonons (those of the pure crystal). This, of 
course, implies that the elements of the resultant scattering matrix, flTTp, are also the sums 
of the corresponding matrix elements for the phonon and impurity scattering. Therefore, 
we write 

p.. = + pW . (21) 

in formulae (0) - ©• 

The electron-phonon contribution p/j^'s have been derived in previous papers, e.g., in 
Ref. H| they are expressed in terms of the Bloch(-Griineisen) functions (H). The appropri- 
ate expressions for the electron-impurity contribution P^ mpS) following from the Yukawa 
potential ([5]) are given in this Section, but first we write down the final expression for the 
electronic thermal conductivity with the influence of the residual electrical resistivity as a 
parameter. For this purpose we have to separate elastic and inelastic contribution to P22 
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and therefore we write explicitly the matrices Uij of the kernels in (10). For the parabolic 
band e(k) = h 2 ~k/2m they can be represented as 

Mll = (^/™) 2 (q-a) 2 , (22) 

W12 = W21 = (ft/™) 2 [(e - + (e' - e)] (q ■ a) 2 , 

M22 = (ft/ m )» [( e - C) 2 + 2(e' - e)(e - C)](q ■ a) 2 + (5' - e) 2 [(q ■ a) 2 + (k • a) 2 ] . 

where (e — e) is the transferred energy during the scattering, q = k' — k is the transferred 
wave vector and terms linear in k and k, which do not contribute to the integrals, have 
been neglected. It is seen that the second component in uyi and the second and third 
one in U22 contribute only when a change of the electron energy in the scattering is taken 
into account. It means that the latter terms describe the effects of inelasticity in the 
scattering. With the definition 

p iph) = p t P h) v (ph) (t) ^ with t = t/Td (23) 

and where Pq^ is a material constant (see Appendix B), we get 

p W = p ( ph ) {KT) 2 [v (i, P H) {t) + V (2, P H) {t)] = + pjW* t (24) 



where 



r 2 



V£> ph \t) = t 5 [\r 7 {t) + ^y 2 Mt) - ^Mt)] , (25) 

and e s is a constant (see Appendix B). 

The explicit form of Py^ (t) describing the electrical resistivity is 

v!t\t) = t s r 5 {t) . (26) 

The two terms in uu, accounting respectively for the elastic and inelastic contributions, 
after integration with the transition probability for the phonon scattering can be repre- 
sented by a single term, 

p (ph) = p( P h) = p (ph) ^ T)V W (t) (27) 



3 These terms describe the basic inelastic contribution. For the phonon scattering the single 
term of uu and the first terms in the case of u%2 and U22 are expressed by the Debye functions 
in case of pure elastic scattering and by the Bloch(-Gruneisen) functions in case of inelastic 
scattering (see Ref.8). The terms proportional to (e' — e) and (e' — e) 2 in u\2 and U22 lead to an 
essential qualitative effect in the thermoelectric power and the thermal conductivity. 
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with 



vi?\t) = p!t\t) = ±i*rt(t) . (28) 

Two material constants influence the magnitude of V^' ph \t). These are the Fermi 
energy, e F , and e s = 2mv 2 , where v s is the sound velocity (see the Table I in Appendix). 
The latter energy (or more precisely e s /4) is a threshold energy for the phonon emission 
in the electron-phonon scattering process (see Appendix J in Ref. |I"8|). 

Assume now that the charged (static) impurity concentration is so low that one can 
consider them to be statistically independent. Then the scattering of the conduction 
electrons by a single impurity contributes solely to the P-ij matrix elements and they 
are linear in their concentrations. In consequence the electrical resistivity is linear in 
the electron concentration and for the degenerate electron gas it is also temperature 
independent. 

The latter statement can be easily confirmed by the calculation of p^ 71 in the 
metallic limit. For further convenience we assume 

Sommerfeld expansion of terms of (^) corresponding to an elastic scattering we have 
p (nnp) = p(irn P ,i) = p ( P h) (^2/3^^2^ Therefore, the final expressions for P n and P 22 
describing the impurity and phonon scattering are 

Pn = P^ h) Pu(t), with P u (t) = Pt ] [r + Pif\t)) 

P22 = Pi ph \k B T) 2 P 22 (t), with P 22 (t) = ^-[r + P<* h) (t)) + P£ ph \t) . (29) 



A. Electrical resistivity and thermal conductivity in metals 

The above expressions for the matrix elements of the P-matrix yield the electrical 
resistivity 

P = Po, Ph [r + P[ p 1 h \t)], (30) 

where po,pft is a constant (see Appendix B), and the electronic thermal conductivity 

t k 2 CoT d 

K — K 0iPh , K ,ph — ~7y • yoL) 

r22{t) o p 0lP h 

Note that the Lorenz number reads 

r-P K -r n2V ^ (30) 

and reduces to the Sommerfeld value C for elastic scattering, when one neglects P22 
in (25). 
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The results (26) describing the electrical resistivity are obviously standard ones and 
do not require any comments. As concerns the electronic thermal conductivity and the 
Lorenz number similar expressions in terms of the Bloch-Griineisen functions as those 
in the present paper were obtained by Wilson! and were numerically examined in Ref. 
|9] manifesting the influence of the impurity scattering. The Wilson expression for the 
thermal resistivity in our notation reads 

T^^ + lrf^^^^ • (33) 



B. Seebeck coefficient of metals. 



Pr 



In order to discuss the thermoelectric power, S, we introduce, as for the discussion 

a relative order of magnitude parameter ip = 
in terms of which for the Yukawa potential (131) we have r(e F ) = <p<P(2k F /\). 



of the thermal conductivity in Ref. 3 I 



(imp) I p (ph) 



It will be also convenient to introduce here the function \P{x) = x 4 /(i + x 2 ) 2 in terms of 
which we express [d$/d£:] eF = (e F )~ 11 P(2k F /X) . In so doing the off-diagonal component 
of the P-matrix responsible for the impurity scattering is linear in k B T/s F and reads 



p(imp) 



(T) 



P n (imp) (fc B T) P$ 



(imp) 



ny(imp) 
1 12 



7T 2 h T 

O Bp 



(34) 



When these expressions are inserted into (§) and linear terms in k B T/e F are only accounted 
for, one obtains the following expressions for the coefficients in (|3"8|) 



S (t) 



7T ks V X2 

3 



( P h) 



(*) 



V 2 2(t) 




7T 3 Vu(t) 
3 V 2 2(t) 



,7T^ 2 y9^(2fc F /A) 



(35) 



Note that the impurity scattering has an influence on the term So (t) only through the 
matrix element V 22 (t) in the denominator. If ip = 0, the above formula corresponds to 
that obtained in Ref. ^] for pure phonon limited thermoelectric power. 

Notice that the pure phonon limited thermoelectric power obtained within the Wilson 
formalismlll do not represent any characteristic features of the experimentally observed S 
behavior in simple monovalent metals. Such results for S obtained in Ref. ||] were already 
discussed there in terms of various approximations. 



IV. TEMPERATURE DEPENDENCE OF THE TRANSPORT COEFFICIENTS 

AND THE LORENZ NUMBER. 

To compute the transport coefficients as a function of t = T /To we assume the value 
of 2k F /X = 0.75 and < <p < 1. The temperature dependence of the reduced electrical 
resistivity following from (26), p/po, P h, which is shown in Fig. 1, and the Seebeck coef- 
ficient in Fig. 2 have been computed for ip = 0.0,0.2,0.5, 1.0,2.0 and the corresponding 
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values of T = <^<2>(0.75) = 0.0, 0.0173, 0.0259, 0.0431, 0.0862, 0.173 respectively. As it is 
seen the impurity scattering slightly lowers the maximum of S following from the transfer 
of the energy between the electron and phonon system!. The values of e s and e F have 
been chosen to be the same as those in Ref. [8|, namely e s = 6.0 K and e F = 1.5 eV. The 
electronic thermal conductivity is considerably affected by the impurity scattering, what 
is seen from the data represented in Fig. 3 and Fig. 4. In Fig.3 the reduced thermal 
conductivity (k/ko,p/i) is shown as a function of t in the low temperature range such as in 
Fig. 4.16 of Ref. ^ where the results of Wilson's considerations!!! were already illustrated. 
In Fig. 4 the results are shown in the same temperature range as the thermoelectric 
power in Fig. 2. We note that he dependence of the thermal conductivity as a function of 
temperature here obtained by the Ziman variational method is qualitatively the same as 
that following from Wilson's considerations!^. Furthermore, at physically realistic tem- 
peratures the values of k obtained by the latter method approach a saturation, which 
is characteristic of monovalent metals. Under the assumption of the same values of the 
parameters describing the impurity scattering the saturation of k following from Wilson 
formula occurs at a so high temperature that it ceases to be meaningful. According to 
ZimanEl, who examined the Wilson formula, physically reasonable results for monovalent 
metals (saturation and lack of minimum at intermediate temperatures) are obtained when 
the Umklapp scattering is taken into account. Experimental results for rare earth inter- 
metallic compounds exhibit the existence of a minimum at temperatures such as those 
seen in Figs. 3 and 4 (see e.g. Fig. 7 in Ref. |2^ and references therein), in which the 
magnetic scattering are also taken into account. 

The reduced Lorenz number, C/Cq, as a function of t is shown in Fig. 5 for the same 
values of the parameters as n/no tP h in Fig. 4. The behavior of this quantity as a function of 
t is qualitatively the same as that following from the Wilson considerations^'!. Note that 
for only phonon scattering £(0) = 0, whereas k(0) tends to infinity at zero temperature. 



V. CONCLUDING REMARKS. 

Using a simple physical electron model, applying the Ziman variational method within 
the Boltzmann equation approach and using the same level of approximations for each 
property we have calculated the phonon-limited electrical resistivity, thermoelectric power 
and the electronic thermal conductivity of a simple metal and studied how the presence 
of charged impurities influences these transport coefficients. The method allows us to 
determine the temperature dependence of the coefficients in a wide temperature interval 
and to express them in terms of simple and standard functions already useful for examining 
experimental data. However we are aware that to take into account a more complicated 
band structure is surely necessary in many cases. In this purpose the Allen methodEl is 
surely more appropriate though the task needs to bear upon numerical technicalities from 
the very beginning. Note that the latter method is strictly applicable to metallic systems, 
when the electron gas is degenerate, while by using the method of this present paper one 
advantage is that one can easily consider semiconductor or semimetallic systems as wellll. 
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A comment is in order on the results of the present paper obtained in the limit of the 
degenerate electron gas. As concerns the electrical resistivity and the electronic thermal 
conductivity (or resistivity) the final results are qualitatively the same as those obtained 
previously in Refs. |l|,|9|,[10|. In the case of the thermoelectric power however we can obtain 
a new evidence as a conclusion, i.e. the maximum at intermediate temperatures is caused 
by the energy transfer between the conduction electrons and phononsa and is lowered as 
a result of the elastic scattering of the electrons on impurities. 

The final formal expression for the electronic thermal resistivity reads as a simple 
expression: this quantity is proportional to the scattering matrix element P22, much like 
the electrical resistivity is proportional to P n . The result is in fact model independent 
and is a consequence of confining the scattering matrix to four elements (see Appendix 
A). 

The range of applicability of the method we have used here above cannot be finally dis- 
cussed without referring to papers using more general and sophisticated methods and/or 
models. In the first place this concerns the papers in which the electron transport was 
studied within the formalism of the Boltzmann equation's In the second place a rela- 
tion to the linear response theoryiBi and the quantum transport equationE^Ej will be 
given. However, by using these methods one cannot obtain simple, compact, and useful 
expressions for the three transport coefficients, moreover being valid in the wide tem- 
perature region experimentally studied. Further complications have also to be expected 
when taking into account the phonon drag@, the virtual recoillllll or the electron-phonon 
r enormalisat ionHS 

Let us start our discussion with a comparison of our model and results with the other 
ones obtained by the Ziman variational method within the formalism of the Boltzmann 
equation. As we stressed at the beginning we confined ourselves to a model of the parabolic 
electron band and the Debye model of the phonon system. This oversimplified model al- 
lowed us to gain a coherent and qualitatively correct description of the three transport 
coefficients: the electrical resistivity, thermoelectric power and the thermal conductivity 
(resistivity) for a simple and standard model of the scattering. Any generalization of 
these results to an arbitrary band structure requires numerical studies at an initial stage. 
A compact formula in this general case of the Ziman variational method has been ob- 
tained only for the electrical resistivity among the three coefficients, a formula which has 
been called "Ziman's resistivity formula"0. It has been used to compute the electrical 
resistivity for a realistic Fermi surface. 

Again a more exact method is that of Allen§ though a numerical one. The electrical 
resistivity is again simplest to compute but the results for the thermal conductivity^ and 
the thermoelectric powerll were also obtained by this method. The essence of the method 
is the separation of the k- variation during the scattering into two types: (i) an "angular 
variation" over the Fermi surface, and (ii) a "radial" variation with e(k) perpendicular to 
the Fermi surface describing the inelastic electron scattering. It is obvious from our results 
for the thermoelectric power! that taking into account inelastic electron scattering, which 
is equivalent to considering a sufficient number of harmonics describing the radial varia- 
tion within the Allen formalism, is of major importance for obtaining a correspondence 
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to experimental results if the scattering of the electrons by excitations such as phonons 
is to be accounted. In Ref. in which the lowest order radial harmonic was considered 
and the analysis was in fact restricted to the Mott formula, a poor correspondence to the 
experimental results for the thermoelectric power was obtained. Inelastic scattering and 
the radial harmonics are also important to describe properly the thermal conductivity. In 
Ref. [14| the radial harmonics were applied to discuss the influence of the electron-electron 
scattering on the thermal conductivity. According to our knowledge no such considera- 
tions were made to analyze the influence of the electron-phonon scattering neither on the 
thermal conductivity nor the thermoelectric power. 

The assumed statistical independence between the impurities as scattering centers 
restrict our considerations to low impurity concentration. This is consistent with the 
experimental fact that the Matthiessen rule is valid below a critical residual resistivity! 
above which the T 5 low temperature dependence (following, for instance, from the Bloch- 
Griineisen dependence of the electrical resistivity) ceases to be valid. One can hardly es- 
timate such a critical concentration without analyzing a specific system and experimental 
results. 

To estimate the magnitude of the transport coefficients out of the reduced values 
illustrated on each figure in Sec. IV we present the order of magnitude of the relevant 
material constants in Appendix B. The magnitude of e s in Table I is the same as in Ref. [| 
and corresponds to the magnitude of the sound velocity and the effective electron mass(es) 
of real materials. 

We start our discussion of more sophisticated methods than the Boltzmann equation 
formalism by noting that in the linear response theory the validity of the Matthiesen 
rule, which was assumed by us in a generalized version, is related to the fact that the 
self-energy of the conduction electrons is the same as the sum of the phonon and impurity 
self-energy^. This holds only in the lowest Born approximation. When higher order 
processes are taken into account, i.e., when interference between the electron-phonon and 
electron-impurity interaction is to be considered, the Keldysh quantum transport equation 
formalismE3c3'L3 is apparently better. 

In Ref. ^ the same model as in the present paper was considered to study the influence 
of the electron-phonon-impurity interference on the thermoelectric power. A considerable 
nonlinearity of S as a function of temperature was found, the stronger the higher impurity 
effective charge. The effect can be hardly compared with that studied in the present paper 
since it has been obtained under an assumption of elasticity for the electron-phonon scat- 
tering whereas the assumption of inelasticity of the electron-phonon scattering is essential! 
for obtaining the maximum of S found in the present paper and Ref. 8. Note that as a 
result of a weak elastic impurity scattering (considered within the Born approximation) 
the bump in S is found to be lowered with the strength of the impurity interaction or the 
impurity concentration. 

The electron-phonon renormalization@@ leads also to a nonlinearity of S as a function 
of T. The effect was studied under the assumption of elastic electron-phonon interaction 
or in the adiabatic approximation. Again there can be no discussion between our results 
and those in Ref. [?| and [30] except for the region of very low temperatures, where the effect 
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of elasticity and renormalization are not effective and where S is a linear function of T 
with the same slope! as in the pure crystalE. 

Though the impurity scattering is known to play an important role in the thermal 
conductivity at least in the low temperature region there is a very limited number of 
papers devoted to this coefficient aspect. The basic difficulties and problems related to 
applying the methods of the linear response formalism to treating the electronic thermal 
conductivity problem in the pure electron-phonon system are described in Ref. How- 
ever, such considerations as those in Ref. [22] have not be generalized for a system with 
impurities. On the other hand the considerations of Refs. [7],[n], devoted to thermoelectric 
effects for a model of Debye phonons with fixed impurities, have not been used by these 
authors in an appropriate way to consider the electronic thermal conductivity. There- 
fore, a complete theory of the electronic thermal conductivity still requires a thorough 
consideration of the mutual role of the electron-phonon-impurity interference, phonon 
renormalization and elasticity of the scattering, much the same as on the thermoelectric 
power. For comparison to data, in not too good metals, the phonon heat contribution 
should be also taken into account for the thermal conductivity. 

A question may arise if our results can be compared to experimental ones in a straight- 
forward, quantitative way. The answer seems that these are the alkali metals, which ex- 
hibit essentially spherical Fermi surface to within a few per cent, that should be first of all 
considered as materials for the comparison. It occurs, however, that the phonon system 
and the electron-phonon interaction in these metals are too complicated to be well ac- 
counted by the applied approximations. It concerns in the first place, (i), the continuous 
medium approximation used to treat the phonon system and the deformation potential 
approximation applied to the electron-phonon interaction (see, e.g., Ref. 12). In the sec- 
ond place, (ii), the Debye approximation is insufficient for the alkali metals. Finally, in 
the description of the electron scattering we also neglect the Umklapp processes, which 
occur important even at very low temperatures0S in the alkali metals. 

The approximations (i) mean that the longitudinal phonons of long wavelength are only 



considered, whereas the transverse phonons^ and those of short wavelength (contributing 
also to the Umklapp process) have to be taken into accountSfi Due to a strong elastic 
anisotropy the Debye approximation (ii) is not appropriate even at low temperatures since 
a constant frequency surface is distorted from the Debye spherical formEae3. Furthermore, 
the strong influence of the electronic system on the behavior of the phonons cause that 
the basic symmetry relations between the elastic constants (Cauchy relations) are no 
longer validfcH and, therefore, even the basic results following from the continuous medium 
approximation of Ref. |T2| may be questioned. All these factors cause that the simple Bloch 
approach have to be used with a great caution in attempts of the quantitative description 
of the electrical and thermal transport in the alkali metals. 

The most early investigations@@ of the alkali metals indicated that one can fix such 
parameters as those of Appendix B to fit the experimental electrical resistivity to the 
Bloch-Griineisen function, provided that the temperature is not too low. Thorough 
analysis0@ based on the measured phonon spectra of the alkali metalsH'S showed later 
that the problem is more complicated and the effects mentioned by us above have to be 
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in fact taken into account. Therefore, even in case of these simplest metals the correct 
quantitative theory based on first principles is in fact impossible and a quantitative anal- 
ysis have to be based on a cross section for the scattering of the conduction electrons 
written in terms of the phonon spectrum resulting from neutron scattering datail'0. Pre- 
cise experimental measurements of the transport coefficients of the alkali metals at low 
temperatures!! support this conjuncture and reveals other physical effects which should 
be taken into account for a proper quantitative analysis in this temperature region. 

As concerns the thermoelectric power our results indicate that other factors than those 
in case of the electrical resistivity are important. The impurity scattering occurred less 
essential than expected. Within the assumed model this is the parameter e s which besides 
e F is solely responsible for the strength of the inelastic scattering of electrons by phonons 
and controls the magnitude and sign of S. The obtained temperature dependence of S is 
the same in character as that obtained experiment ally@ for all the the alkali metals except 
Li (which exhibits a peculiar behaviori but the Fermi surface of which comes closest to 
the Brillouin zone boundaries). In search for a quantitative theory of S(T) of a specific 
metal one should also estimate the order of the contribution of the phonon dragp, that 
following from the electron-phonon renormalization0 and the electron-phonon-impurity 
interference. The described above factors which influence the electrical resistivity in the 
alkali metals may occur of secondary importance. 

The effects of the impurity scattering manifest in most pronounced way in the behavior 
of the thermal conductivity but due to the validity the Wiedemann-Franz law in the limit 
of low temperatures the effects of the impurity scattering can be easily extracted out of 
our results at low temperatures for a comparison with experimental ones23. As concerns 
the comparison of the values of appropriate coefficients following from the theory and 
experimental data for the alkali metals the same problems as in case of the studies of the 
electrical resistivity may arise. The Umklapp process are expected to be less dominant 
at low temperatures than in case of the electrical resistivity but the electron-electron 
interactions may occur essential for the thermal conductivity since both the experimental 
and theoretical investigations indicate^ that electron-electron scattering in the alkali 
metals impedes heat flow much more effectively than charge flow. 
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APPENDIX A 



The transport coefficients considered in the paper follow from the Onsager expressions 
for the electrical J and thermal U currents: 

J = aE + LVT, 
U=ME + iVVT, (36) 

which in the "experimental" form are following 

E = p J + 5 VT , 

U = 7T J - k VT (37) 

with 7r = — ST and 

p=l/a, S = -pL , k = MpL - N . (38) 
Within the Ziman variational method the coefficients of (|36|) are 

& ^ ] JiQijJj i L — ^ ] JiQijUj , 

i,j=l i,j=l 

M=J2 U iQn J j > N = -h E ^^'^ » ( 39 ) 

i,j=l i,j=l 

where J, and Z7j (i,j = 1,2, • • •) are respectively the electron and thermal trial currents 
defined by (p!2[)-(p^) and Qij is the inverse matrix to the scattering matrix P^-, i.e. Qij = 
(P _1 )ij. In the approximation consisting in abbreviating the combinations in (|39|) to two 
components (z,j = 1,2) we obtain (7)-(9) from fl3"gp. For a degenerate electron gas, i.e. 
for sufficiently large z = e F /k B T, the final expressions for the transport coefficients are 
(14) - (16). Notice that the thermal resistivity in particular is proportional to P22 • 
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APPENDIX B 



The functions describing the trial currents, as written below (15)-(16), are 
L (z) = (3/2) F l/2 (z) , 

Li(z) = -(3/2)zF 1/2 {z) + (5/2)F 3/2 (z) , (40) 
L 2 (z) = (3/2)z 2 F l/2 (z) - 5zF 3/2 (z) + (7/2)F 5/2 (z) . 

For completeness we also write down their asymptotic forms 

L (z) - z 3 ' 2 , Lx(z) = (vr 2 /2) z 1 ' 2 , L 2 (z) = (tt 2 /3) z 3 ' 2 . (41) 

The coefficients of the trial currents are Jq = e Uq with 

MW 2 ^) 2 ( v - a ) 

Uq — { when the factor is (42) 

1 /2mA 5/2/ 7^ 



i.e. the upper case corresponding to the factor (v • a), and the lower one to (k ■ a) in the 
trial function (|]). 

The coefficient which determines the dimension and the order of the magnitude of the 
phonon scattering matrix depends on the effective electron mass, m, the volume of the 
primitive cell, V, its total mass, M, the sound velocity v s , the phonon wave vector cutoff, 
qn, and the deformation potential constant, E x . It reads 

■ 487r3Mi, 3 (^) 2 (m) 2 ^ ( V " a ) 

P = { 2 when the factor is (43) 

487t3Md 3 (#) 2 ^ ( k-a )' 

The coefficient po,p/i = -Pq /(e 2 [/g£ F 3 ) in fl3U| ) can be represented by 

Po, P k = 48n(ri 2 (^)Q)VE 2 ~p 0!ph (44) 
where m /(k 2 e 2 ) = 9.113 • 10 3 [eV" 2 nm~ 3 ], 

Po, ph = {^-f (vsgo)- 1 = (^-f h(k B T D y\ (45) 

where is the free electron mass and p = m/m Q . For the value of Tjj = 2007^ assumed 
for our representative data we have h(k B T[))~ l = 3.437 ■ 10~ 2 Q cm. In Table I we show 
some representative data of po^h for e F = 1.5 eV, T D = 200 K, p = 1 and the values of e s 
given in the table. They correspond to those in Ref. |||, where numerical relations between 
e s and other constants describing the phonon system are presented. 

The screened Coulomb potential scattering contribution to the transport is described 
in terms of the dimensionless functions $ and ^ defined by (^) and fl54"|). The coefficient 
p(imp) q £ scattering matrix elements is 
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(im) \skn, mp Z^CmB 2 (v-a) 
p(tmp) _ i when the factor is (46) 

i^^V(f) 2 (k-a). 

where nj mp is the impurity concentration (their number per unit volume). The residual 
resistivity p<* mp ) = P 1 ^ mp - ) /(e 2 f/Qe P 3 ) can be represented by 

p (mp) = ^^p imp <P(2k F /X) (47) 
n e p 

where n e is the number of the conduction electrons per unit volume, with 

Pimp = 3 ; (48) 

m VF 

vf is the Fermi velocity, and 2k F /\ = 1 as we mentioned in Sec. II. For the data 
represented in Table I (free electron gas, e F = 1.5 eV) we have pi mp = 0.075 f2 cm. 

In order to relate the magnitude of r in (25)-(26) with the number of impurities in 
the sample we represent it as <p 

r = P - = ^ <?(2A; F /A) = ^<P(2k F /\) (49) 

Po, P h Po, P h n e n n e n 

and shall estimate (p. It means that (p (31) is ip — (Z 2 n imp /n e p)<p. For the mass of 
the unit cell M = 10 6 m , its volume V = 0.1 nm 3 and E\ = (2/3)e F = 1 eV we have 
p 0ph = 0.1 po, P /i- This is an underestimated value in many cases, thus we may assume 
Po,ph — a Po,ph with 0.1 < a < 1. For e F = 1.5 eV, To = 200 K p, = 1 the value of p imp is: 
p imp = 0.1 Qcm. The latter value along with those ones for p ,p/i following from the values 
of po, P h presented the Table I yields an estimate of (p. 

Finally, for information as to the magnitude of the thermal conductivity represented 
in dimensionless units in the figures we give the magnitude of 

7T 2 CqT D 

K , P h = — (50) 

<J P0,ph 

for the Debye temperature To = 200 K in the last column of Table I. 
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TABLES 



TABLE I. The values of the parameters allowing to estimate the magnitude of the electrical 
resistivity and electronic thermal conductivity for a free electron effective mass, when the Fermi 
energy is e F = 1.5 eV and the Debye temperature Tr> = 200 K. 



no. 


e s [K] 


po, p h [fi cm] 




Ko, P h [W/mK] 


(i) 


1.50 


1.93 -10~ 3 


10° - 10 1 


2.54 -10~ 5 


(ii) 


3.00 


2.39 -10~ 4 


10 1 - 10 2 


2.04 -10~ 4 


(iii) 


6.00 


2.99 -10~ 5 


10 2 - 10 3 


1.64 -10" 3 


(iv) 


12.00 


3.79 -10- 6 


10 3 - 10 4 


1.23 -10~ 2 


(v) 


24.0 


4.58 -10- 7 


10 4 - 10 5 


1.07 -10- 1 


(vi) 


48.0 


5.76 -10- 8 


10 5 - 10 6 


8.48 -10- 1 
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Figure Captions 



Figure 1 Temperature dependence of p/po iP h resulting from Eq.(26) for values of 
the parameters 2k F /\ = 0.75 and the relative scattering strength ip = 0.0,0.2,0.5,1.0. 
They correspond to the values of F = ^(0.75) which are 0.0, 0.0173, 0.0431, 0.0862, 
respectively. 

Figure 2 — Temperature dependence of S/Sk, e for values of the parameters as in Fig.l 
with e ¥ = 1.5 eV, e s = 6 K, and the effective electron mass equal to that of a free electron. 
S k , e = k B /e = 86fi V/K. 

Figure 3 — Temperature dependence of the reduced thermal conductivity below T D , 
for 2k F /\ = 0.75 and ip = 0.0, 0.2, 0.5, 1.0, 2.0; the corresponding values of F are 0.0, 
0.0173, 0.0431, 0.0862, 0.173, respectively. For F = 0, the value of re is infinite. Note 
that for the magnitude of r as large as 2.0 there is a plateau of re instead of the maximum 
at T = 0.1Tb. 

Figure 4 — Temperature dependence of the reduced thermal conductivity in the same 
temperature range as that displayed for the thermoelectric power in Fig. 2. The value of 
2k F /\ is the same as in Fig.3, i.e. 0.75, and tp = 0.0, 0.2, 0.5, 1.0, 2.0. 

Figure 5 — Temperature dependence of the reduced Lorenz number in the same tem- 
perature range as that of the thermoelectric power in Fig.2. The value of 2k F /\ is 0.75 
as in Fig.3 , and tp = 0.0, 0.2, 0.5, 1.0, 2.0. 
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